TITLE OF THE INVENTION 

METHOD FOR ITERATIVE DETERMINATION OF DISTANCE BETWEEN 
RECEIVING STATION AND TRANSMITTING STATION AND ALSO 
CALCULATING UNIT AND COMPUTER SOFTWARE PRODUCT 

CROSS REFERENCE TO RELATED APPLICATIONS 

[0001] This application is based on and hereby claims priority to Gemian Application No. 
103 16 290.9 filed on 9 April 2003. the contents of which are hereby incorporated by reference. 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

[0002] The invention discloses a method for the iterative determination of the distance 
between a receiving station and a transmitting station and also a corresponding calculating unit, 
a corresponding computer software product and a corresponding data carrier. 

2. Description of the Related Art 

[0003] The GPS method (GPS: Global Positioning System) for detemnining the position of a 
station is based on the principle of measuring signal transfer times between satellites having 
known positions and a GPS receiver that is installed In a station. The transfer time for a signal is 
proportional to the distance of the respective satellite from the station, that is, the distance 
between a satellite and the station can be converted, using the propagation velocity of the 
corresponding signal, into the respective signal transfer time. Since both the position of the 
station and the extent by which the clock at the station deviates from the atomic clocks on the 
synchronized satellites are unknown, the position of the station is detenmined using the 
distances from at least four satellites. 

[0004] The signals which the satellites transmit are satellite-specific code sequences which 
are modulated to a carrier frequency. The basic measurement of the station to determine its 
position consists in measuring for each satellite a phase offset between the respective code 
received and a copy of the code, which is produced in the receiver The respective phase offset 
(chip code offset) is a measure of the signal transfer time and therefore of the distance between 
the station and the conresponding satellite. A code available for civil use currently consists of 
1023 chips and is repeated with a cycle of one millisecond. In the one-millisecond duration of a 
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code cycle, a satellite signal which is propagated in a vacuum at the speed of light travels about 
300 km. The distance between the satellite and the ground station is typically several thousand 
kilometers, however. During the time which the signal requires to cover the distance between 
the satellite and the station, however, the code is therefore repeated several times. In other 
words, this means that, during the signal transfer time, the satellite uses the code N times to 
modulate the continuously transmitted signal and that the station receives the code N times in 
the same period of time. 

[0005] A precise determination of the phase offset is an essential component of a 
determination of the distance that the station is from the respective satellite and therefore forms 
the basis of the determination of the station's position. Since the phase offset Is, however, 
measured using only one code cycle, the total transfer time that the signal requires to cover the 
distance between the satellite and the station cannot be determined from the phase offset 
alone. The signal transfer time is made up of a whole number of repeats of the code and of a 
modulo ensuing from the determination of the phase offset. The phase offset is calculated 
directly by the station, whereas the integral multiples of the code are calculated by an iterative 
method in the course of the determination of the station's position. The computing power that 
such an iterative method needs is directly dependent on the number of iterative steps required • 
to determine the integral multiples in a satellite's code. 

SUMMARY OF THE INVENTION 

[0006] The invention therefore addresses the problem of providing an improved method of 
determining the distance between a satellite and a station. 

[0007] In the method according to the invention for the iterative determination of the distance 
between a receiving station and a transmitting station, information is received by a calculating 
unit via a signal from the transmitting station that is received by the receiving station, the signal 
being provided with a code. A number of integral multiples of the code come into the receiving 
station during the signal transfer time that the signal requires to travel between the transmitting 
station and the receiving station. During the first iteration the number of integral multiples in the 
code is calculated and the calculated number is used in the second iteration at least without a 
fresh calculation being earned out. As a result of the fact that in the second iteration at least, a 
fresh calculation of the number of integral multiples is not made, the computing power required 
can be reduced. 
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[0008] It is therefore particularly advantageous if the calculated number of integral multiples 
is used in all the iterations. 

[0009] It is advantageous that a rough estimation of the position of the receiving station is 
used to calculate the number of integral multiples of the code. The rough estimation of the 
position of the receiving station, together with an estimation of the position of the transmitting 
station, allows a sufficiently precise estimation of the distance between the receiving station and 
the transmitting station to be made. From the distance calculated in the above way it is possible 
to calculate the exact number of integral multiples if the position of the transmitting station is 
known almost exactly (in satellites, for example, from their known stable orbits) and if the 
position of the receiving station is known with an accuracy greater than the half-wavelength of 
the code. Here, the wavelength of the code is obtained by multiplying the duration of the code 
by the propagation velocity of the signal, that is. in vacuum, by the speed of light. The exact 
number of integral multiples is consequently detemriined in the very first iteration of the method 
and is no longer necessary in all subsequent iterations. When multiplied by the wavelength of 
the code, the number of integral multiples corresponds to a minimum distance between the 
transmitting station and receiving station. The above minimum distance is one of the addends 
used to calculate the actual distance in the first and at least, in the second iteration. 

[0010] It is practical if the receiving station is located in a radio cell of a radio communication 
system and if the rough estimation of the position of the receiving station is based on a radio cell 
identifier assigned to the receiving station. As a result of such an estimation of the position, the 
position of the receiving station Is known with an accuracy greater than the half-wavelength of 
the code. The exact number of integral multiples can be calculated in the very first iteration 
therefore, as described in the previous paragraph. 

[0011] The number of iterations required to detemnine the distance between the receiving 
station and the transmitting station can be reduced, that is, the computing power required can 
be reduced, if a value other than zero is used in the very first iteration for the signal's signal 
transfer time. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0012] These and other objects and advantages of the present invention will become more 
apparent and more readily appreciated from the following description of the exemplary 
embodiments, taken in conjunction with the accompanying drawings of which: 
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Figure 1 is a block diagram of an exemplary embodiment of a radio system according to 
the invention which shows the transmission of signals by four satellites to a mobile station and 

Figure 2 is an illustrative diagram of a signal transmission between a satellite and the 
mobile station shown in Figure 1 . 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

[0013] Reference will now be made in detail to the preferred embodiments of the present 
invention, examples of which are illustrated in the accompanying drawings, wherein like 
reference characters refer to like elements throughout. 

[0014] A receiving station is any station that is able to receive signals. In the text that 
follows, a mobile station is considered to be a receiving station. A mobile station is, for 
example, a mobile phone or also a device the location whereof can be moved, for the 
transmission of image and/or sound data, for sending faxes, for sending short messages via a 
Short Message Service (SMS), for sending emails and for Internet access. It is therefore a 
general transmitting and/or receiving unit in a radio communications system. 

[0015] A transmitting station is understood to be a station that is able to transmit a signal 
provided with a code, i.e. a modulated signal. A transmitting station may be assigned to a radio 
communications system or be an external station that exchanges data with the radio 
communications system or is used by the radio communications system. In particular, a 
transmitting station is a satellite in a GPS system. 

[0016] Signals that are provided with a code or modulated are both GPS signals and any 
other signals used in relation to a detemiinatlon of distance or determination of position. 

[0017] The invention can be used advantageously in any radio communications systems. 
Radio communications systems are understood to be any systems wherein a data transmission 
ensues between stations via a radio interface. The data transmission can ensue both bi- 
directionally and uni-directionally. Radio communications systems are in particular any mobile 
radio systems, for example those based on the GSM (Global System for Mobile 
Communication) or the UMTS (Universal Mobile Telecommunication System) Standard. Future 
mobile radio systems, for example fourth generation systems, are to be understood as being 
radio communications systems. 
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[0018] In the text that follows, the invention is described using as an example a mobile radio 
system based on the UMTS Standard, without however being restricted thereto. 

[0019] Figure 1 shows a diagram of a signal transmission by four satellites SAT1 , SAT2. 
SAT3, SAT4 to a mobile station MS. which is located in a radio cell Z. A cell identifier ZID of the 
radio cell Z is assigned to the mobile station MS. Via the cell identifier ZID, a rough 
geographical position, based on the known geographical position of the radio cell Z, is assigned 
to the mobile station. The con^esponding geographical position is typically the center of the 
radio cell Z, the position of a base station BS that supplies the radio cell Z, or even the weighted 
mean of a plurality of radio cells, if the mobile station MS is at the same time connected to a 
plurality of radio cells. 

[0020] The four satellites SAT1 . SAT2, SAT3. SAT4 transmit to the mobile station MS the 
signals SI, S2. S3. S4 respectively, which are modulated with a satellite-specific code C1, C2, 
C3, C4. The signals SI, S2, S3, S4 have a respective signal transfer time t1. t2, t3. t4 and are 
used to detemriine the distance between the mobile station MS and the con-esponding satellite 
and to determine the position of the mobile station MS. Position data for the satellites SAT1 . 
SAT2, SAT3, SAT4 is known to the mobile station MS. The mobile station knows the orbits of 
the satellites SAT1 , SAT2, SAT3, SAT4 from the corresponding Almanac AL and ephemeris 
information EP, which is transmitted thereto from a radio access network UTRAN by a base 
station BS. The mobile station MS evaluates the signals S1, S2, S3. S4 received by the 
satellites SAT1. SAT2. SAT3 and SAT4 and. for each signal SI. S2, S3, S4, it determines a 
phase offset of the respective code C1 . C2, C3. C4 received with respect to the respective 
reference codes produced in the mobile station MS. The reference codes concur with the codes 
C1 . C2. C3, C4 for the respective satellites SAT1 , SAT2. SATS. SAT4 and are produced by the 
mobile station MS in such as way that, if the mobile station MS and the satellites SAT1 . SAT2. 
SAT3, SAT4 were to be in perfect synchronization, there would be no phase shift if the 
corresponding satellite SAT1 . SAT2. SATS, SAT4 and the mobile station MS were to be located 
in the same place. The phase offset that is actually measured is obtained from the transfer time 
t1. t2, tS, t4 of the corresponding signal S1, S2, S3, S4 of the satellites SAT1, SAT2. SATS, 
SAT4 and from an imperfect synchronization of the atomic clocks on satellites SAT1 . SAT2, 
SATS, SAT4 with one another and with a clock on the mobile station MS, which achieves time 
measurement for the mobile station MS. 
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[0021] Signals S1 , S2, S3, S4 of the satellites SAT1 , SAT2, SAT3, SAT4 contain respective 
so-called ephemeris infomnation, which records inter alia a deviation of the respective satellite 
SAT1, SAT2, SATS, SAT4 from its stable orbit and contains infonmation on the en^or in the clocks 
on satellites SAT1, SAT2, SATS, SAT4. The mobile station MS transmits information I to the 
base station BS of the radio access network UTRAN, which infonnation contains in particular 
the phase shifts measured for codes C1 , C2, CS, C4. The cunrent ephemeris infonmation is 
known in the radio access network via a control GPS receiver, for example. Of course, the 
infonmation I can also contain cun-ent ephemeris infonmation for satellites SAT1 , SAT2, SATS, 
SAT4, should the data not exist in the UTRAN radio access network. 

[0022] The base station BS is connected to a radio network controller RNC, wherein a 
calculating unit BE is located. In the calculating unit BE, an iterative method of detenmining the 
distance between the mobile station MS and the satellites Is carried out. From the satellite 
distances, the position of the mobile station MS is further detemriined. The method uses the 
Almanac information and the ephemeris infonmation from satellites SAT1, SAT2. SATS, SAT4 
and also the phase shifts of the codes C1 , C2, CS, C4 from the satellites SAT1 , SAT2, SATS, 
SAT4, the shifts having been measured by the mobile station MS. The radio network controller 
RNC, which is connected to a core network CN of the radio communications system, is able to 
pass on to the core network CN the position detenmined by the calculating unit BE. From the 
core network CN, the position can then be transmitted to any further communications systems. 

[0023] Of course, the calculating unit BE can also be ananged in the mobile station MS or in 
the core network CN. 

[0024] Figure 2 shows a diagram of the transmission of a signal S1 from a satellite SAT1 to a 
mobile station MS. The signal S1 is modulated with a code C1. which has a duration of tC1 
seconds and which Is repeated accordingly every tC1 seconds. The duration tC1 of the code 
C1 is distinctly shorter than the total time t1 that the signal S1 requires to be transmitted over 
the distance D1 between the satellite SAT1 and the mobile station MS. 

[0025] The signal transfer time t1 can be expressed by the duration tC1 of a code cycle. 
Hence 

t1 = tC1 * (n, + ei). 
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[0026] In the above formula, ni is a whole number ^ 0 and conresponds to the number of 
integral multiples of the code C1 produced during the signal transfer time t1 . The in-ational 
number 81 is between 0 and 1 and con^esponds to fragments of a code C1 . 

[0027] For the distance D1 between the satellite SAT1 and the mobile station MS, the 
following formula applies: 

D1 = t1 * V, 

V being the propagation velocity of the signal S1 . 

[0028] The phase offset between the code C1 received and the reference code produced in 
the mobile station MS is the in'ational number 81. It is detemnined directly by the mobile station 
MS. The integral multiples ni of the code C1 and also the clock error on the satellite SAT1, and 
the clock enror on the mobile station MS are determined by the iterative method in the 
calculating unit BE. 

[0029] The invention has the advantage that the first iteration only is used to detennine the 
integral multiples ni. Thus a saving is achieved compared to existing methods, with respect to 
the calculations and computing power required. By a cell identifier ZID of the radio cell Z in the 
mobile station MS, the geographical position of the mobile station MS is known with greater 
accuracy than is the half-wavelength A of the code C1 prior to the detemriination of position 
using GPS. Here, the wavelength A of the code C1 is equal to the product of the propagation 
velocity v of the signal S1 and the duration tC1 of the code C1. Where the duration tC1 of the 
code C1 is for example one millisecond, the half-wavelength A of the code C1 is about 150 km. 
Since typical cell sizes in mobile radio systems are in the order of size of 20 km, a knowledge of 
the cell identifier ZID means that the position of the mobile station MS can already be 
determined more readily than the half-wavelength A even at the beginning of the detemiination 
of position. The number of integral multiples n^ of the code C1 can therefore be calculated 
correctly in the very first iteration according to the invention. A fresh calculation of the number 
of integral multiples ni is no longer necessary in the following iterations. 

[0030] Individual calculations of the method are desaibed below with reference to Figure 1 
and Figure 2: 
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[0031] A system of coordinates for calculating the position is, for example, selected in such a 
way that the origin of the system of coordinates is in the center of the earth. Such a system of 
coordinates is known in English as Earth-centered-Earth-fixed (ECEF). In the above exemplary 
embodiment, the codes C1. C2, C3. C4 that are used all have a duration tC1 of one millisecond 
and each consists of 1023 chips. The invention is not of course restricted thereto but can also 
be applied forthwith to other durations of the code sequences C1 , C2, C3, C4 and accordingly to 
different numbers of chips. The invention can likewise be used where the codes C1 , C2, C3, C4 
do not all have the same duration or the same number of chips. ^ 

[0032] The phase offset between the codes C1 , C2, C3, C4 received and the respective 
reference codes produced in the mobile station MS is determined by the mobile station MS and 
expressed in the fomi of whole chips Xwj and fragments Xf,i of chips received. For further use in 
the method, the calculating unit BE converts the phase offset into a distance equivalent Cj. 

Ci = A(Xw,i+2-'°X,.,)/1023 

[0033] The subscript i stands for the respective satellite SAT1 , SAT2, SATS, SAT4 and A 
stands for the wavelength of the respective codes C1 , C2, C3, C4 which, in the above 
exemplary embodiment, is the same for all the codes C1, C2, C3, C4. 

[0034] The method for determining the distance now calculates the position of the respective 
satellite SAT1 , SAT2, SAT3, SAT4 at the time of transmission of the corresponding signal S1 , 
S2, S3, S4. Only the respective time of reception is known, however, with the result that in the 
first iteration of the method, an estimation of the time of transmission is then made for the 
respective signal S1, S2, S3, S4. As a starting value for the iterations, a time of transmission U 
is selected, the time being for example ti° = 67.48 milliseconds prior to the time of reception 
te • ts" = te - ti° applies. Here the superscript stands for the number of iterations, beginning with 
zero In the first iteration. The above value is given by the fact that the satellites SAT1 , SAT2, 
SAT3, SAT4 are at a distance of at least 20,230 km away when in a vertical position above the 
mobile station MS located on the earth's surface. In conjunction with the propagation velocity v 
of the signals SI , S2, S3, S4. a minimum signal transfer time of 67.48 milliseconds can thus be 
calculated and be preset as a starting value for the method in the very first iteration. Using the 
value for the time of transmission that has been estimated in this way, a first estimated position 
of the respective satellite S° is calculated. 
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[0035] Of course, starting values other than tj° = 67.48 milliseconds can also be selected. 
However, any corresponding starting value tj^is always based on a realistic assumption relating 
to the minimum distance between satellites SAT1 . SAT2, SATS. SAT4 and the mobile station 
MS. Consequently, different starting values tj^can also be selected for the different satellites 
SAT1. SAT2, SAT3. SAT4. 

[0036] Furthermore, to determine the distance, the method requires the en-or Atj*^ of the 
respective satellite clock. The superscript k is the iteration index and begins with k=0 in the first 
iteration. The method is able to calculate the satellite clock enror Atj*" with the aid of a clock 
correction vector at, which is contained in the ephemeris information. 

At(' = ao' + aV (te - ti^ - toe') + 32' (t - t-' - toef 

ao', B^ and 32 are the components of the clock correction vector aj and toe' is the reference time 
for the ephemeris information. 

[0037] The method uses the error At|^ to correct the distance equivalent Cj of the phase shift. 
i'^ = Ci + cAti'^ 

[0038] In a subsequent calculation of the method, the integral multiples nj of the codes C1 , 
C2, C3, C4, produced during the transfer time t of the respective signal S1, S2, S3. S4. are 
calculated. To do this, the method uses a rough estimation of the position of the mobile station 
MS. It uses the cell identifier ZID of the radio cell Z, the identifier being assigned to the mobile 
station MS, in order to estimate accordingly the geographical position of the mobile station 
MS. Assuming that the above position is more accurate than the half-wavelength A of the 
codes C1 , C2, C3, C4. in other words more accurate in the above exemplary embodiment than 
about 150 km. the exact value of the integral variables nj in the first iteration, I.e. for k=0, is: 



ni = int{(||R°-S|° 



II- I* 



i°)/A} 



int(x) being the next whole number to the number x. with int(x) ^ x. 



[0039] In subsequent iterations, i.e. for k 



> 0. ni is not re-calculated. 
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[0040] The clock error on the mobile station MS is the same for all the satellites SAT1 , SAT2, 
SATS, SAT4 and is therefore calculated using data from one satellite only. The following 
formula generally applies: 

C= i' + niA-||R°-Si''|| 

[0041] Using the integral variables n-, of the codes CI . C2, C3. C4, the distance equivalent 

i\ of the phase offset and the clock error tb" of the mobile station MS, corrected by the satellite 
clock en-or Atj\ a distance pi"" from the mobile station MS, known as the pseudo range, is 
calculated for each satellite SAT1, SAT2, SATS, SAT4. 

Pi''= i'^ + niA+ctb" 

[0042] Using as a basis the distances p,'' of the satellites SAT1 , SAT2. SATS, SAT4, a new 
estimated value R""*^ is determined for the position of the mobile station MS. By the new 
estimated value and of the position Si*" of the respective satellite that has been estimated for the 
above iteration, a new estimated value for the signal transfer time ti*"*^ is calculated between the 
respective satellite SAT1, SAT2, SATS, SAT4 and the mobile station MS and used in the 
following iteration. The following fomiula applies: 

t.k*l = ||Rk*-'_S>||/c 

[0043] The method is considered to be converged if the following condition has been met: 
II R^*^ - r"" II 5 

[0044] In the above case no further iteration is canried out. The distances pi"" between the 
mobile station and the satellites are fixed and the current estimated position R"*^ is established 
as an actual position and converted into an output format that gives the geographical latitude, 
longitude, and altitude of the mobile station MS. The uncertainty of the position is at 
Au = II R''*' - R"* ||. 

[0045] If the condition || R"*^ - II > ^2 is met, however, then the method interrupts any 
further calculations and issues an error report. 
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[0046] The invention has been described in detail with particular reference to preferred 
embodiments thereof and examples, but it will be understood that variations and modifications 
can be effected within the spirit and scope of the invention. 
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